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^) , The properties of spin polarized neutron matter are studied both at zero and finite temperature 

(N 

Q ' using the Dl and the DIP parameterizations of the Gogny interaction. The results show two 

(D ' 
Q ■ different behaviors: whereas the DIP force exhibits a ferromagnetic transition at a density of 

Pc ~ 1-31 fm~^ whose onset increases with temperature, no sign of such a transition is found for 

^ \ Dl at any density and temperature, in agreement with recent microscopic calculations. 
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The possible existence of a phase transition of neutron matter to a ferromagnetic state 
has motivated many investigations of the Equation of State (EoS) of spin-polarized neutron 
matter. Besides the interest that such a transition could have in the context of neutron stars 

3een addressed in the framework of 



[l|, this problem has gained interest in itself and has 
very different theoretical approaches [2, Isl, \4, M, la, lZI, Isl, M, llO, luj- Whereas some of these 
calculations, like for instance those based on Skyrme-like interactions, predict a transition 
at densities in the range (1 — 4)po (with po = 0.16 fm~^ the saturation density of symmet- 
ric nuclear matter), others, like recent Monte Carlo [7| or Brueckner-Hartree-Fock (BHF) 
calculations [9|, [lO|| using modern two- and three-body realistic interactions, exclude such 



a transition at least up to densities around five times po- In spite of this discrepancy, it 
is interesting to study how temperature influences the ferromagnetic transition [l2|, [l3|. In 



Ref. [l^ an analysis of the temperature effects in the framework of a Hartree-Fock calcula- 
tion with Skyrme interactions was performed. In the present paper, we would like to study 
the influence of the flnite-range terms of the interaction on the ferromagnetic transition. 
To this end, we will consider the Gogny interaction. This is an effective nucleon-nucleon 
force with both zero- and finite-range terms and a simple spin-isospin structure. In addition 
to the original Dl parameterization [1^, several other parameterizations of this force have 
been obtained. The DIS force, for instance, was introduced to improve the pairing prop- 



erties and surface effects of finite nuclei [15|, while more recently the DIP 



16l | interaction 



has been introduced with the aim of reproducing the EoS of pure neutron matter given by 



a variational microscopic calculation with realistic interactions |17|. 

The properties of nuclear matter deduced from Gogny interactions have already been 
treated in the literature ll8|. Indeed, several instabilities produced by these forces at zero 

. . n . . . .. 

temperature have been studied in previous works [19|. The isospin instability, for instance, 
is a common feature to all the existing Gogny parameterizations as well as of most Skyrme 
forces. This instability is signaled by the fact that, above a certain critical density, the energy 
per particle of nuclear matter becomes more repulsive than that of neutron matter. For the 
DIP force, this instability takes place at pi ~ 7po, while for Dl it occurs at pi ~ 3po. 
Furthermore, Dl and DIS also exhibit a spinodal instability in neutron matter, i.e., the 
energy per particle does not increase monotonically with density. Instead, it reaches a 
maximum value at a critical density (which for Dl and DIS is around ps ~ 4po) and then 
decreases from that density on. For the discussion of our results, we will choose the DIP 



and Dl parameterizations of the Gogny NN force, since they represent two quahtatively 
different behaviors regarding ferromagnetic instabihties. 

In general, any Gogny interaction can be casted in the following form: 

2 2 

Vnn{^ = J2 ^~^ (^* + B,P^- H, Pr - M, P, P,) + to (l + xo P.)p% 5{f) (1) 

with r the distance between two nucleons. The spin-isospin structure of the force is given by 
the spin (isospin) exchange operators P^ {Pt)- Notice that the usual spin-orbit term of the 
Gogny interaction has been omitted, because it does not give any contribution in infinite 
matter. The Gogny force includes a sum of two gaussian-shaped terms which mimic the 
finite-range effects of a realistic interaction in the medium. Usually, it also contains one 
density-dependent zero-range term, even though in the case of the DIP parameterization 
two of these terms were used to make the fitting procedure more flexible. 

In the following, we will consider spin-polarized neutron matter, characterized by the spin 
polarization parameter A = (p| — Pi)/{p^ + Pi), where p|(|) is the density corresponding 
to neutrons having spin up (down) respect to a given direction, and by the total density, 
P = Pt + Pi- The energy per particle in the Hartree-Fock approximation is given by 



1 v^ h^k^ 
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+ — 2^ {ki<7i,k2cr2 I Vnn I ki(Ti,k2cr2)A x n„^{ki,T) n„^{k2,T) , (2) 
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where n„{k,T) are the momentum distributions: 

of each spin component a. At each temperature T and density p„^ a self-consistent proce- 
dure has to be performed to compute the chemical potential of each species, p^, from the 
normalization condition: 

Pa = Y,n^{k,T). (4) 

k 

Once the momentum distribution is determined, the entropy per particle of the system is 
given by: 



s(PT,Pi,^) = -^ {n^{k,T)\nn^{k,T) 

P cr,k 

+ [l-n.(A;,T)]ln[l-n.(A;,T)]}. (5) 



From the internal energy and the entropy one can readily compute the free energy per 

particle, f{pi^,pi,T) = e{p^,pi,T) — Ts{p-i, pi,T), and, from this, the inverse magnetic 

susceptibility is given by: 

1 1 f Pi^f\ 
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with /i, the magnetic moment of the neutron. 

An important quantity for our analysis is the single-particle (sp) energy: 

5e h^e 
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where Ua{k) is the sp potential. Note that in contrast to the sp spectrum e^ik) appearing 
in Eq. ([3]), this sp energy contains the rearrangement effects. The momentum dependence 
of the sp potentials can be characterized by the effective mass m* (A;): 
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Let us start our analysis with the properties of the sp potential, U{k). Fig. [T] reports the 
sp potential obtained from DIP for neutrons with spin up at saturation density in a non- 
polarized system (A = 0) and in a fully polarized one (A = 1) for T=0 (left panels) and 40 
MeV (right panels). The splitting due to the different spin polarizations can be understood 
in terms of both the dependence of the sp potentials in phase space and the dependence on 
spin of the effective NN interaction 13|. Moreover, the increase of temperature makes the 
sp potentials slightly less attractive in both cases due to the dependence on the transferred 
momentum of the two-body matrix element which, at finite temperature, is explored in a 
wider range of momenta. In spite of the fact that the effective mass depends on momentum, 
for low enough values of k the sp potential can be casted in a quadratic form: 

1 dU„{k) 



Ui'\k) = U^{0) + 



2k dk 
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in which the derivative 



dU^jk) 
dk 



and the ^-independent term f/cr(0) are determined from the full 



sp spectrum at the corresponding temperature. Fig. [T]indeed illustrates that, for the density 
and polarization under study, this quadratic approximation (dotted lines) is very close to 
the full fc-dependent U(j{k) sp potential (full lines), at least for momenta up to kp. The 
approximation of Eq. ((91) is devised to reproduce the slope of the spectrum at kp together 



with its value at k 



0. The choice of the Fermi momentum kp is related to the fact that the 



most relevant changes in the occupation of the sp states, when changing p and T, essentially 
involve the neighboring region of the Fermi surface. The quadratic approximation of Eq. ([9]) 
leads to an approximated form for the sp energies in terms of the effective mass at the Fermi 
surface and a k-independent term of Ua yk): 

Finally, let us note that a similar picture is obtained for the other Gogny parameterizations. 
In Fig. [2] we report the ratio XFree/x (where XFree is the magnetic susceptibility of the 
free Fermi sea at the corresponding temperature) as a function of density. Results for the 
DIP (Dl) forces are shown on the left (right) panels of the figure at several temperatures. 
A different qualitative behavior for the two parametrizations is clearly observed. On the one 
hand, the DIP parameterization leads, similarly to what was found in the case of the Skyrme 



interaction 



12l |. to a ferromagnetic phase transition, signaled by a vanishing ratio XFree/ X- 
The critical density of this transition at T=0 MeV is pc = 1.31 fm~'^, a much larger value 
than the ones obtained with Skyrme interactions, all of which were systematically below 1 



fm~^ (see Tab. II of Ref . [l2|)- Moreover, and even though it is difficult to distinguish this 
in the figure, the critical density slightly increases with temperature. This is in agreement 
with the intuitive idea that thermal disorder increases the onset density of ferromagnetism, 
but it is opposite to what was found with Skyrme forces. On the other hand, no trace of a 
ferromagnetic transition is seen for Dl at any density or temperature. This behavior is very 



similar to what was found in a BHF calculation with realistic interactions [13|. In addition, 
the temperature dependence (the higher the temperature, the lower the ratio becomes) is 
also in agreement with such microscopic calculations. 

In this context, it is interesting to study the behavior of the entropy. To this end, plots of 
the entropy per particle as a function of spin polarization at a fixed value of density (p = Po) 
and several temperatures are shown in Fig. [3l We consider the entropy given by Eq. ([5]) 
computed with both i) the full momentum-dependent spectrum ecr{k) (symbols), and ii) the 
quadratic approximation to ea^{k) of Eq. ((91) (lines). In the two cases the exact chemical 
potential determined from the normalization of the momentum distribution, Eq. (Ill), are 
used. The agreement between the values of the entropy coming from the exact and the 
approximated sp energies is rather satisfactory, the most significant discrepancies not being 
larger than a 3%. Both the Dl and the DIP parametrizations give rise to similar results: the 



entropy per particle is symmetric with respect to the non-polarized state and shows a maxi- 
mum at zero polarization. In addition, the values of the entropy increase with temperature, 
which is again an indication that the entropy per particle behaves as naivel y ex pected. This 
so-considered "natural" behavior was also found in the BHF analysis of Ref. 13|. In contrast, 
for Skyrme forces the entropy per particle of the polarized phase is seen to be higher than 
the non-polarized one for a certain density on [l2|. This defines a kind of "critical" density, 
which is smaller than po for most Skyrme parametrizations. Such a non-intuitive behavior 
of the entropy as a function of the polarization can be related to the dependence of the 
entropy on the effective mass, and a condition for the effective masses can then be derived: 

^*iP^^ = ^) < 2^/3 fll) 

m*(p,A = 0) ' ^ > 

if the quantity s(p, A = 1) — s(p, A = 0) has to be always negative. Most of the Skyrme 



forces analyzed in Ref. [12|| violate this criterion and thus lead to an "anomalous" temperature 
dependence for the onset density of ferromagnetism. 

Notice, however, that for Skyrme forces the effective mass is momentum and temperature 
independent. This is actually not the case of neutron matter described by means of Gogny 
forces: the effective masses do depend on both momentum and temperature. Nevertheless, as 
it has been previously discussed, the sp spectrum and the entropy per particle are correctly 
described by a quadratical momentum dependence, with an effective mass calculated at 
k = kp for each temperature. Within this approximation, one can proof that the criterion of 
Eq. (ITTl) is still valid at each temperature. We have checked that both the Dl and the DIP 
forces fulfill the criterion of Eq. (ITTl) in a vast region of the (T, p) parameter space, which 
includes both the classical and degenerate limits. Thus, as it has been previously observed, 
we do not expect that the ferromagnetic transition (only present for DIP) has an anomalous 
thermal behavior. 

One of the questions that needs for a closer insight is why neutron matter described with 
Gogny forces does not present a ferromagnetic instability, or why it is only present at very 
high densities. Intuitively, one expects that the instability will be related to the zero-range 
term. This can be understood just by taking into account that the pure zero-range term 
only involves S-wave contributions and, therefore, due to Pauli principle, it is not possible 
to have a couple of neutrons with the same spin interacting through the contact term. 
Therefore there is no contribution to the total energy per particle of fully polarized matter 

6 



from a zero-range term, whereas for non-polarized neutron matter this term is both strongly 
density-dependent and repulsive, thus contributing to the spin instability. In contrast, the 
behavior of finite-range terms, both direct and exchange, has the opposite dependence with 
polarization: for a given density, the higher the polarization, the higher these contributions 
become. In other words, if only finite-range terms were considered, the non-polarized system 
would be energetically favored. In addition, the density dependence of these contributions 
is softer than that of the zero-range term. Therefore, the competition between zero and 
finite-range effects is resolved in favor of the first one at high values of the density. From 
this, it is also obvious to understand why no transition is found for the Dl and the DIS 
parametrizations. In these two cases, the fitting parameters are such that the zero-range term 
does not contribute for any polarization, because Xq = 1. The zero-range term is therefore 
not active in neutron matter, whereas the direct and exchange parts of the finite-range 
terms behave as mentioned, energetically favoring the non-polarized case at all densities. In 
contrast, for the DIP force xq ^ 1, and the zero-range term is active. As a consequence, the 
instability arises once we come across a high enough critical density. 

In this work, we have studied the properties of polarized neutron matter with neutrons 
interacting through Gogny forces, both at zero and finite temperature. The results show two 
different qualitative behaviors for the two parametrizations under study. On the one hand, 
the DIP exhibits a ferromagnetic transition at what can be considered a very high density, 
Pc ~ 1.31 fm~^. On the other hand, no sign of a ferromagnetic transition is found for the Dl 
parameterization at any density or temperature. The cause of these two different behaviors 
is related to the zero-range term, which contributes for DIP, but not for Dl. Finally, we 
have checked that the temperature dependence of the ferromagnetic transition with DIP is 
the one expected by intuition, in agreement with the fact that the criterion of Eq. (fTTl) for 
the effective masses is respected in a wide range of densities and temperatures. 
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FIG. 1: Single-particle potential V(}z) for up neutrons in non-polarized (top panels) and totally 
polarized (bottom panels) neutron matter at p = 0.16 fm^^ for T = (left panels) and T = 40 
MeV (right panels) within the DIP parameterization. C/(/c) and its quadratic approximation V^'^' 



are displayed in full and dashed lines respectively. The arrows denote the value of kp. 
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FIG. 2: Ratio between the inverse magnetic susceptibility of interacting neutron matter and that 
of the corresponding free Fermi sea as a function of density for several temperatures. 
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FIG. 3: Entropy per particle as a function of the spin polarization at p = 0.16 fm~'^ and several 
temperatures. The exact Hartree-Fock values are depicted with circles, diamonds and triangles, 
whereas the lines illustrate the results obtained from a quadratic approximation of the single- 
particle spectrum (Eq. (fTOl) ). 
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